Extended JC-Dicke model for a two-component atomic BEC inside a cavity 
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We consider a trapped two-component atomic Bose-Einstein condensate (BEC), where each atom 
with three energy-levels is coupled to an optical cavity field and an external classical optical field 
as well as a microwave field to form the so-called A- type configuration. After adiabatically elim- 
inating the atomic excited state, an extended JC-Dicke model is derived under the rotating-wave 
approximation. The scaled ground-state energy and the phase diagram of this model Hamiltonian 
are investigated in the framework of mean-field approach. A new phase transition is revealed when 
the amplitude of microwave field changes its sign. 
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I. INTRODUCTION 

The so-called Dicke model [l| describes the interaction 
of a large number of two- level systems (e.g., atoms) with 
a single optical mode. Since the effective light-matter 
coupling strength is dependent on the number of atoms 
N (oc \/~N), a sufficiently large N would lead to a classical 
phase transition 0, |3| (at finite temperatures) from the 
normal state, which corresponds to the atomic ground 
state associated with the vacuum state of the optical 
mode, to the superradiant-phase state, where the phe- 
nomenon of superradiance occurs with the finite scaled 
mean numbers of both photons and excited-state atoms. 
Recently, exploration of quantum phase transitions in the 
Dicke model at zero temperature has attracted significant 
attentions [3-[2^|. The drastic change across the critical 
point due to a qualitative change of the ground state of 
the Dicke model has been investigated, in the frameworks 
of the scaled ground state energy, macroscopic atomic 
excited-state population, quantum entanglement, Berry 
phase, quantum chaos and so on. 

Apart from the standard Dicke model, several more 
generalized Dicke models have also been proposed and 
studied 0-i, In some of them [H EH, the 

free atoms in the standard Dicke models are replaced 
by the atomic Bose-Einstein condensate (BEC). A BEC 
describes a collective quantum state of a large number of 
atoms and may be used to generate a macroscopic quan- 
tum object with a longer lifetime compared to the free 
atoms. It has attracted much interest to combine the 
BEC with optical cavity in the strong coupling regime 
[23l - l26j . On the other hand, comparing with the case 
of free atoms in the standard Dicke model, there is an 
additional atom- atom collision interaction in the BEC, 
which could lead to the change of critical phenomenon 
[l5l ]. More recently, under some approximations, Chen et 
al. [l(| proposed an extended Dicke model based on a 
two-component atomic BEC in an optical cavity with 
the atoms coupled to both the quantized optical cavity 
filed and an external classical driving optical field, where 
an interesting phase diagram covering phenomena from 



quantum optics to atomic BEC is addressed. 

In this paper, we consider a trapped two-component 
atomic BEC where the two condensated states and an 
ancillary excited state of each atom form a three-level A- 
type configuration via the couplings to the optical cavity 
and an externally controlled classical optical field as well 
as an external microwave field. Such a three-level config- 
uration can be reduced to a two-level configuration (i.e. 
the two condensated states) after adiabatically eliminat- 
ing the atomic excited state under the large sigle-photon 
detuning condition. Under the so-called rotating-wave 
approximation (RWA), we derive an effective, extended 
JC-type [29| Dicke model for such a two-component BEC. 
The phase diagram for the derived extended JC-Dicke 
model is also investigated in detail. 



II. MODEL AND HAMILTONIAN 

The setup of the JC-Dicke model under consideration 
is depicted in Fig. [TJ An optically-trapped Rb atomic 
BEC under the two-mode approximation with the atomic 
states 5 2 Si/ 2 \F = l,m/ = —I) (ground state |1)) and 
|F = l,m/ = 0) (metastable state |2)) is placed in a 
single-mode quantized optical cavity. These two states 
of the atomic BEC are coupled via both an external mi- 
crowave field and a two-photon process 27j mediated by 
an ancillary excited state |3) (from h 2 Pj,/2), where the 
single-photon transition between |3) and the ground state 
1 1) (or the metastable state |2)) is coupled to the quan- 
tized optical cavity field (or the external classical opti- 
cal field). Here we assume that both the corresponding 
single-photon dctunings are large and the correspond- 
ing two-photon detuning is very small. For such a case 
of two-photon Raman process 



281 ] . the ancillary excited 



state can be adiabatically eliminated and the effective 
Hamiltonian of the two-component BEC reads (H = 1 
hereafter) 



H, 



ph 



Hat-ph + Hat-cl + Hat-c 



(l) 



2 



with 

Hat = V\c\ C l + ( V 2 + ^2 - LOi)c\c 2 , 

H ph = ua'a, 
Hat-op = A e ge tWcl *4cia + h.c, 

Hat-rnw = fle'^^dd + k.C, 



(2a) 
(2b) 
(2c) 
(2d) 



Hat-at = ycJcjciCl + -^-c\c\c2C2 + 7?12Ci C l4 C 2 ( 2 e) 

denoting the energies of the free atoms in the BEC, the 
free cavity field, the reduced effective interaction be- 
tween the BEC with the optical fields, the interaction 
of the BEC with the microwave field, and the atom- 
atom collision interaction of the BEC, respectively. Here, 
c i.2 (c\ 2) are the annihilation (creation) bosonic opera- 
tors for |1) and |2), respectively; cjj (i = f,2, 3) is the 
corresponding internal level-energy for atomic state 
Vi = Jd 3 r^*(r)[-V 2 /2m + y(r)]^(r) (I = 1,2) is the 
trapped frequency for the states |1) and |2) with V(r) be- 
ing the trapped potential, m the atomic mass, and (j)i(r) 
the corresponding condensate wavefunction; a (at) is the 
annihilation (creation) operator of the cavity mode with 
the frequency to; A c ff = <7i3^23/A is the reduced effective 
coupling strength for two-photon Raman process, where 
<?13 (^23) is the corresponding coupling strength between 
the BEC and the quantized cavity field (classical optical 
field), A is the large single-photon detuning: A = W3— 
(i>2 + UJ2) ~ w ci ^ {313^23} with ui c \ the frequency of 
the classical optical field; fi is the corresponding cou- 
pling strength between the BEC and the microwave field 
(with the frequency w mw ); 77; = (Anpi/m) J d 3 r \4>i(r)\ 4 
and 7/12 = (47rpi 2 /m) / d 3 r |0*(r)</> 2 (r)| 2 with p t and p 12 
(— P21) the intraspecies and the interspecies s-wave scat- 
tering lengths, respectively. It is remarked that the RWA 
has been used for all optical/microwave fields coupling to 
the atomic BEC. 

By using the Schwinger relations 



J+ = cJ,ci, 



■h = 



c\c2 



J_ = c\c 2 , 
cjci 



(3a) 
(3b) 



which fulfill 



[J+,J_] =2J Z , [J Z ,J±] = ±J±, (4) 
the Hamiltonian (TTJ) can be written as 

H cS = (data + ujqJ z + jzJl 



me 



where 



+ -^=e iWclt a) J+ + h.c] + const, (5) 



(6) 
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FIG. 1: (Color online) (a) Schematic diagram of the exper- 
imental setup for a trapped BEC of 87 Rb atoms, (b) The 
internal atomic level configuration of BEC: The |1 )«->-|2) tran- 
sition of the BEC atoms is coupled to the external microwave 
field with the coupling strength Q; The |l}-0-|3) and |2)«-»|3) 
transitions are coupled largely-detuned to the optical cavity 
and the external classical optical field with the correspond- 
ing coupling strengths (713 and respectively, (c) When 
the single-photon detuning A is large (the corresponding two- 
photon detuning is assumed to be small) , a reduced two-level 
configuration can be obtained after eliminating the ancillary 
excited state |3). A c ff is the effective coupling strength for the 
two-photon Raman process. 



is the number of the atoms, 

N-l , 

w = v-2 + Lo 2 - v\ - u\ H — (772 - 771) , (7a) 



v = ( — 2 

A = Acffv^, 
and the constant term 



(7b) 
(7c) 



N 7/ 2 77l 

const = — [(v 2 + u 2 - cJi - — + 72 N) + (ui- — + 7i7V)J 

can be neglected in the following consideration. 

For D, 7^ 0, to eliminate the time-dependence in Hamil- 
tonian we perform a unitary transformation U = 
exp[— iu} mw J z t — i(id mw + Lo c \)cvat\ and obtain an effec- 
tive Hamiltonian 

H = (d a ata + oj b J z + jrJz 



+ (^=aJ + + OJ+ + h.c), 
WN 



N = c\c 2 + c\ci 



(8) 



where u> a = (lj — id mw — uj c \) and oj = (uj — cj mw ) are 
the effective frequencies in the rotating frame, in which 
H is independent of time. To the best of our knowledge, 
Hamiltonian (|8|) appears to be a new one in literatures, 
and thus we call it as the extended JC-Dickc model. 



III. MEAN-FIELD GROUND STATE ENERGY 

We now look into the ground-state properties of Hamil- 
tonian ((HJ) and the corresponding quantum phases as well 
as their transitions. We here consider the case of posi- 
tive u> a , where the stable ground state is anticipated for 



3 



HamiltonianJHJ) . By using the Holstein-Primakoff trans- 
formation [30j ] 



bWN-tfb, J- = VN - tfbb, 



2 



(9) 



with [b,^] = 1, the effective Hamiltonian reads 

H = U a a?a + U b (tfb - y) + ^(6+6 - y ) 2 



[(Aa + Q^)6t A /i_^ 



(10) 



Similar to that used in the standard Dicke model 
we here introduce the displacements for the two shifting 
boson operators as c) = a) — \J~N a* and S = b^ + ^/N (3* 
with the complex displacement parameters a and (3 de- 
scribing the scaled collective behaviors of both the atoms 
and the photons 0, H, [IH, Il6| , In fact, the current 
method of introducing the displacements is equivalent 
to the mean field approach. In this framework, it is clear 
that < |/3| < 1. 

After expanding the terms in the square in (|10p . the 
scaled Hamiltonian can be written up to the order of TV -1 
as 



H/N = H + N~ x / 2 Hi + N~ 1 H 2 , 



(11) 



where 



Ho = tu a a*a + u; b (j3* /3 - i) + r)(/3*/3 - \? 



[(Aa + n) py/l-F0 + cc] 



(12) 



denotes the scaled constant energy, and Hi% denote the 
linear and bilinear terms, respectively. It is noted that 
Hq,i,2 are independent of the number of atoms N. 

The scaled ground-state energy is just given by the 
scaled constant energy in the Hamiltonian 



Eg (";P) = Jq - H 0, 



(13) 



where the displacements a and (3 should be determined 
from the equilibrium condition 



d[E g (a,P)/N]/da* = 0, 
d[E g {a,0)/N]/dp* =0. 

After some derivations, we find that a is given by 

A* 



a = — /?v/l - /?*/?. 



(14a) 
(14b) 



(15) 



and f3 satisfies 



n^l- I3*p = /3[uj b + w (2/3*13 - 1) 
Vt(3* 



+ ( 



2^1- (3* [3 



+ c.c.)], 



(16) 



where w = r\ + |A| /u) a . 

From the above equation (|16p . it is obvious that J7//3 
should be real since all of the parameters, except for the 
coupling strengths A and fi, are real. Here, we assume 
the case of real O for simplicity [3l|. That means, (3 
should also be real: — 1 < /3 < 1 and satisfy 



= u b Py/l - I 32 + O (2/3 2 -l)+wp (2/3 2 - 1) Vl-/? 2 - 

(17) 

For a general real fi, the scaled ground state energy in 
Eq. (fl"3"j) is given by the displacement (3 (by using Eq. (fT5|) 
to eliminate the displacement a) as 



E?(J3) = M/3 2 -I)- 2n/3VT^+w(f3 2 - l -f. (18) 



The displacement (3 is the non-trivial real solution for 
equation (|17[) . In general, there are more than one real 
solutions for Eq. (fT7|) . only the one that leads to the 
minimal scaled ground state energy should be chosen. 

We would like to remark that the displacements deter- 
mined from the equilibrium equations could just make 
the linear term Hi be 0. At the same time, the bilin- 
ear term Hi makes no contribution to the scaled ground 
state energy. Therefore, the exact forms of H\ 2 are not 
needed in our analysis, and thus only the constant term 
Hq determines fully the scaled ground state energy of the 
current system. 



E„ /w 




FIG. 2: (Color online) (a) The atomic displacement j3 for 
the ground state, (b) the square of atomic displacement /3 2 
for the ground state (or the scaled magnetization plus 1/2: 
m+1/2), and (c) the scaled ground state energy versus Wb 
for different coupling strength Q. The energy and frequencies 
are in units of w (w > 0). 
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IV. NUMERICAL RESULTS AND ANALYSIS 

In this section, we focus on the numerical calculations 
of the scaled ground state energy, e.g., E^ (/3) in Eq. (fTHj) . 
and the corresponding displacement (3 that makes Eg (f3) 
minimal. The minimum E*(fi) (as well as the corre- 
sponding f3) is determined by the three parameters: wt>, 
f2, and w. 

Figure [2] plots (3, f3 2 , and Eg '(/?) (corresponding to the 
ground state of the JC Dicke model) as a function of 
ojb for several values of £1. All the energies/frequencies 
in Fig. [5] are in units of the positive w. We mentioned 
"positive w" here since w may be positive or negative, 
while the negative w will lead to different results. Seen 
from Fig. [21 a second-order normal-superradiant phase 
transition at ujb /w = 1 should be expected in the absence 
of microwave field (ft = 0), as in the standard Dickc 
model 3 ; while in the presence of microwave field (O ^ 
0), the normal phase and the corresponding transition 
disappear. 
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FIG. 3: (Color online) The atomic displacement f3 for the 
ground state versus the coupling strength Q for different u)b- 
The frequencies are in units of w (w > 0). 




FIG. 4: (Color online) (a) The atomic displacement f3 for 
the ground state and (b) the scaled ground state energy Eg 
versus Q and w. The frequencies/energies are in units of LOb 
(here we keep > 0). 

The displacement (3 for the ground state is plotted as 
a function of f2 (in units of positive w) for several val- 
ues of u>b in Fig. [3] For a positive/negative f2, the cor- 
responding (3 is also positive/negative. At the point of 
O — > 0, the ground-state f3 may have a jump. In order 
to look the possible transition phenomenon at the point 
of il — > 0, we plot the ground-state (3 and Eg ((3) as the 



function of the parameters £1 and w (in units of positive 
tOb) in the 3-dimensional (3D) Fig. [4] In experiments, 
the parameters f2 and w are controllable, e.g., the former 
one can be easily controlled by changing the strength of 
the microwave field, the latter one can be controlled by 
adjusting the atom-atom interaction interactions by the 
magnetic- field or/and optical Feshbach resonance tech- 
niques [32M341 ] . From Fig. HJa), it is clear that the dis- 
placement f3 has a jump when w > Ub at the point of 
Q —> 0. Notably, the scaled ground-state energy is al- 
ways continuous at the jump point, but its first deriva- 
tive with respect to O does not, which implies a new kind 
of the first order phase transition at the parameter point 
whenever Q changes its sign. Making the replacement: 
O — > —51, we can find the corresponding replacements 
(3^-/3 and Egjfi) E$(-0) = E* \(3) according to 
Eqs. dTTJ) and (ffS]), just as seen in Fig. HI 




FIG. 5: (Color online) The atomic scaled magnetization M 
versus Q and w. The parameters are in units of \ub\ for posi- 
tive (a) and negative (b) u^. 



(a) (co„>0) (b) K<0) 




-2-1012 -2-101 



w/co f w/\m b \ 

FIG. 6: (Color online) The Q-w phase diagram for (a) the 
positive LOf, and (b) negative u^. 

We may define the scaled 'magnetization' as M = 
(J z )/N = f3 2 — 1/2 for the ground state. A positive 
magnetization M means the atomic inversion population: 
more atoms stay in the upper state |2) than those in |1), 
while a negative M means the opposite case. Especially, 
M = —1/2 corresponds to the normal phase: all atoms 
stay in the lower states. Figure [5] plots the magnetization 
against the parameters f2 and w, where Figs.^a) and (b) 
correspond respectively to the positive and negative Wf,. 
From Fig. [SJ we may define phases Pi, 2, 3, 4, as denoted in 
Fig. [5] and Table |U The line Lq of normal phase and L12 
of the superradiant phase, appearing only in the absence 
of microwave field, are the board lines of the phases 
in Fig. I^a); so are L' and L34 the board lines of 
in Fig. HJb). Table U shows the differences of relevant 
quantities for all phases in Fig. |S] [HJ. 
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Wc now address the points/line of quantum phase tran- 
sition in the phase diagram. Point A in Fig. |5{a) de- 
notes the "well-known normal-supcrradiant phase tran- 
sition (along the transverse arrow) in the standard JC 
Dicke model, while point D represents a similar one (cor- 
responding to point A) for the case with the atomic near- 
inversion population when Wf, is negative. The line L12 
corresponds to the phase transition line segment from Pi 
to P2, which is of the first-order, as indicated above. In- 
tuitively, the phases Pi and P 2 may be viewed as the 
para- and dia- "magnetic" phases, because ^ > and 
< for Pi and P2 , respectively. Correspondingly, the line 
segment of L34 is the first-order phase transition line that 
separates P3 and P4, which may be viewed respectively 
as the dia- and para- "magnetic" phases. 



TABLE I: Relevant quantities for different phases Pi, 2, 3, 4 in 
the ground state. 



quantity 


PI 


P2 


P3 


P4 


OJb 


> 


> 


< 


< 


ft 


> 


< 


> 


< 


P 
M 


(-l.o) 


(-^,0) 

(-1,0) 


(^,1) 
(0,1) 


(-!,-#) 

(o,|) 


BE" 

an 


< 


> 


< 


> 


dM 
dw 


> 


> 


< 


< 


dM 

an 


> 


< 


< 


> 



It is worth pointing that we have neglected the anti- 
resonant terms and the corresponding A 2 terms (A is the 
vector potential of the optical field) in the original Hamil- 
tonian. In the standard Dicke model, it was pointed 
out that the anti-resonant terms would also bring an un- 
neglectful influence near the critical point and lead to the 
modification of the result of the quantum phase transi- 
tion [3, [36[ . It was also pointed out that the quantum 
phase transition in the standard Dicke model happens in 
the effective ultra-strong matter-light coupling regime, 



where the A 2 term could also become very strong and 
may not be omitted. Notably, if the effect of the A 2 
term were not neglected, the quan tum phase transition 
would be impossible to happen [37| in the standard Dicke 
model. In order to obtain the quantum phase transition 
in the Dicke model, it was proposed in Ref. [l2[ to get 
an effective Dicke model which does not include the A 2 
term. In the current scheme, similar to Ref. (l2| . we have 
obtained an effective extended JC-Dicke model that does 
not have the anti-resonant terms and the A 2 terms. In 
our original Hamiltonian, it is assumed that the matter- 
light couplings are much smaller than the corresponding 
atomic transition/optical carrier frequencies, then it is 
safe to neglect the anti-resonant terms and the A 2 terms. 
After the unitary transformation, the time-independent 
effective Hamiltonian may lead to the quantum phase 
transition when the matter-light coupling is comparable 
to the effective carrier frequencies. 

V. CONCLUSION 

In conclusion, we have derived an extended JC-Dicke 
model for a two-component BEC coupled to the quan- 
tized optical cavity and the external classical optical field 
as well as a microwave field. The scaled ground-state en- 
ergy and the phase diagram of this model Hamiltonian 
have been investigated in the framework of mean- field ap- 
proach. A new first-order phase transition has also been 
revealed when the amplitude of micromavc field changes 
its sign. 
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